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ABSTRACT 

We study the expected distribution of massive black hole (MBH) spins and its evolution with cosmic 
time in the context of hierarchical galaxy formation theories. Our model uses Monte Carlo realizations 
of the merger hierarchy in a ACDM cosmology, coupled to semi-analytical recipes, to follow the merger 
history of dark matter halos, the dynamics of the MBHs they host, and their growth via gas accretion 
and binary coalescences. The coalescence of comparable mass holes increases the spin of MBHs, while 
the capture of smaller companions in randomly-oriented orbits acts to spin holes down. We find that, 
given the distribution of MBH binary mass ratios in hierarchical models, binary coalescences alone do 
not lead to a systematic spin-up or spin-down of MBHs with time: the spin distribution retains memory 
of its initial conditions. By contrast, because of the alignment of a MBH with the angular momentum of 
the outer accretion disk, gas accretion tends to spin holes up even if the direction of the spin axis varies 
in time. In our models, accretion dominates over black hole captures and efficiently spins holes up. The 
spin distribution is heavily skewed towards fast-rotating Kerr holes, is already in place at early epochs, 
and does not change much below redshift 5. If accretion is via a thin disk, about 70% of all MBHs 
are maximally rotating and have radiative efficiencies approaching 30% (assuming a "standard" spin- 
efficiency conversion). Even in the conservative case where accretion is via a geometrically thick disk, 
about 80% of all MBHs have spin parameters o/wbh > 0.8 and accretion efficiencies >12%. Rapidly 
spinning holes with high radiative efficiencies may satisfy constraints based on comparing the local MBH 
mass density with the mass density inferred from luminous quasars (Soltan's argument). Since most 
holes rotate rapidly at all epochs, our results suggest that spin is not a necessary and sufficient condition 
for producing a radio-loud quasar. 

Subject headings: black hole physics - cosmology: theory - galaxies: active - galaxies: nuclei - quasars: 
general 



I. INTRODUCTION 

Dynamical evidence indicates that massive black holes 
(MBHs) reside at the centers of most nearby galaxies 
(Richstone et al. 1998). The tight correlation observed 
between their masses and the stellar velocity dispersion or 
mass of the host bulge (Gebhardt et al. 2000; Ferrarese & 
Merritt 2000; Haring & Rix 2004) suggests a close rela- 
tionship between the growth of black holes and spheroids 
in galaxy halos. It is not yet understood which physical 
processes established this correlation and how it is main- 
tained through cosmic time with such a small dispersion. 

Besides their masses, tobh, astrophysical black holes are 
completely characterized by their spins, S — aGmBn/c, 
< s/mBH < 1. The spin of a MBH is expected to 
have a significant effect on its observational manifesta- 
tion. For example, the spin determines the efficiency of 
converting accreted mass into radiation and has implica- 
tions for the direction of jets in active nuclei (Rees 1984). 
The electromagnetic braking of a rapidly spinning black 
hole may extract rotational energy (Blandford & Znajek 
1977), convert it into directed Poynting flux and electron- 
positron pairs, and power some radio galaxies and gamma- 
ray bursts. The orientation of the spin is thought to deter- 
mine the innermost flow pattern of gas accreting onto Kerr 
holes (Bardeen & Petterson 1975). The coalescence of two 



spinning black holes in a radio galaxy may cause a sudden 
reorientation of the jet direction, perhaps leading to the 
so-called "winged" or "X-type" radio sources (Merritt & 
Ekers 2002). 

The spin of a MBH is determined by the competition be- 
tween a number of physical processes. Black holes forming 
from the gravitational collapse of very massive stars en- 
dowed with rotation will in general be born with non-zero 
spin (e.g. Fryer, Woosley, & Heger 2002). An initially 
non-rotating hole that increases its mass by (say) 50% by 
swallowing material from an accretion disk may be spun 
up to o/tobh = 0.84 (Bardeen 1970). While the coales- 
cence of two non-spinning black holes of comparable mass 
will immediately drive the spin parameter of the merged 
hole to o/tobh 0.8 (e.g. Gammic, Shapiro, & McKin- 
ney 2004) , the capture of smaller companions in randomly- 
oriented orbits may spin down a Kerr hole instead (Hughes 
& Blandford 2003, hereafter HB). 

In this paper we make a first attempt at estimating the 
distribution of MBH spins and its evolution with cosmic 
time in the context of hierarchical structure formation the- 
ories. In the next section we review the assembly of MBHs 
in popular cold dark matter (CDM) cosmogonies. We then 
address the evolution of MBH spin via coalescences (§3) 
and gas accretion (§4). We compute the spin distribution 
of MBHs in §5, and summarize our results and discuss 
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their implications in §6. 



2. ASSEMBLY AND GROWTH OF MBHs 

The main features of a plausible scenario for the hi- 
erarchical assembly, growth, and dynamics of MBHs in 
a ACDM cosmology have been discussed in Volonteri, 
Haardt, & Madau (2003), and Volonteri, Madau, & Haardt 
(2003). Dark matter halos and their associated galax- 
ies undergo many mergers as mass is assembled from 
high redshift to the present. The halo merger history is 
tracked backwards in time with a Monte Carlo algorithm 
based on the extended Press-Schechter formalism. "Seed" 
holes form with intermediate masses in the rare high-cr 
peaks ( "minihalos" ) collapsing at z — 20 — 25 (Madau 
& Rees 2001), and their growth and dynamical evolu- 
tion is followed in detail with a semi-analytical technique. 
The merging - driven by dynamical friction against the 
dark matter - of two comparable-mass halo+MBH sys- 
tems ( "major mergers" ) drags in the satellite hole towards 
the center of the more massive progenitor, leading to the 
formation of a bound MBH binary with separation of ~ 
pc. The long dynamical frictional timescales leave many 
MBHs "wandering" in galaxy halos after a merger between 
halos of large mass ratio ( "minor mergers" ) . 

A two-component model for galaxy halos is adopted, 
where the dark matter is distributed according to an NFW 
profile (Navarro, Frenk, & White 1997), while the initial 
central stellar distribution is a singular isothermal sphere 
with stellar velocity dispersion, c*. The latter is related to 
the halo circular velocity following Ferrarese (2002). The 
present-day mass density of nuclear MBHs accumulates 
via gas accretion: in every major merger, the heavier hole 
accretes at the Eddington rate a gas rest mass Amo- This 
leads to a change in the total mass-energy of the hole given 
by 

Am = 3.6 x 10 6 Mq V c % , (1) 

where V Ci i5o is the circular velocity of the merged system in 
units of 150 kms -1 (Volonteri, Madau, & Haardt 2003). 
The scaling with the fifth power of the circular velocity 
of the host halo and the normalization are fixed in order 
to reproduce the observed local itibh — c* relation. The 
quantities Am and Amo are related by Am = (1 — e) Amo, 
where e is the mass-to-energy conversion efficiency, equal 
for thin-disk accretion to the binding energy per unit mass 
of a particle in the last stable circular orbit. 4 This model 
was shown in Volonteri, Haardt, & Madau (2003) to re- 
produce well the observed luminosity function of optically- 
selected quasars in the redshift range 1 < z < 5. 

2.1. MBH binaries 

The semimajor axis of a bound MBH binary continues 
to shrink owing to dynamical friction from distant stars 
acting on each hole individually, until the pair becomes 
"hard" when the separation falls below ah — Gm2(4<7*) 
(Quinlan 1996), where m2(< mi) is the mass of the lighter 



hole. After this stage the MBH pair hardens via three- 
body interactions, i.e., by capturing and ejecting at much 
higher velocities the stars passing within a distance of or- 
der the binary separation (Bcgclman, Blandford, & Rees 
1980). Our scheme assumes that the "bottleneck" stages 
of binary shrinking occur for separations smaller than 
ah', during a galactic merger, after a dynamical friction 
timescale, we place the MBH pair at ah and let it evolve. 
The hardening of the binary modifies the stellar density 
profile, removing mass interior to the binary orbit, de- 
pleting the galaxy core of stars, and slowing down further 
hardening. The above scheme can be modified by rare 
triple black hole interactions, when another major merger 
takes place before the pre-existing binary has had time 
to coalesce (Saslaw, Valtonen, & Aarseth 1974). In this 
case there is a net energy exchange between the binary 
and the third incoming black hole, resulting in the ejec- 
tion of the lighter hole and the recoil of the binary. The 
binary also becomes more tightly bound. Our scheme im- 
plicitly neglects the depopulation of the "loss cone," i.e. 
it assumes a large supply of low-angular momentum stars 
(as one might expect in significantly flattened or triaxial 
galaxies, Yu 2002). It also neglects the role of gaseous - 
rather than stellar dynamical - processes in driving the 
evolution of a MBH binary (e.g. Escala et al. 2004; Ar- 
mitage & Natarajan 2002). 

If stellar dynamical and/or gas processes drive the bi- 
nary sufficiently close (< 0.01 pc), gravitational radiation 
will eventually dominate angular momentum and energy 
losses and cause the two MBHs to coalesce. 

3. SPIN CHANGE BY BINARY COALESCENCES 

HB have recently addressed how binary coalescence 
changes the spin of the remnant hole. We follow here their 
treatment. The problem is simpler for small mass ratios, 
q = mijmx -C 1, as the binary is then well described by a 
test particle, the orbit of which shrinks due to gravitational 
wave emission to the last stable orbit (LSO): the smaller 
hole then plunges into the more massive one adding the 
LSO angular momentum to its spin (since a black hole 
spin scales with its mass squared, the small hole's spin 
can be neglected). While the test particle description is 
not accurate for q > 0.5, binary coalescences with mass 
ratios larger than 0.5 take place in only 10% of our sim- 
ulations at 1 < z < 5, and are even more rare at z < 1 
(see Figure 1). Moreover, the test particle description is 
qualitatively correct in that q > 0.5 coalescences lead to 
rapidly spinning holes most of the time (see Figs. 2 and 
3). 

We adopt below geometrized units and set c = G = 1; 
hatted quantities have been made dimensionless by di- 
viding out powers of mass, e.g. a\ = a\jm\. Be- 
fore coalescence, the larger hole has mass mi and spin 
|S| = aim\. Kerr geodesies orbits are specified by choos- 
ing their energy E, angular momentum parallel to the 
spin, L z , and "Carter constant" Q. Orbits with Q ^ are 
inclined at angle 9 with respect to the equatorial plane, 



cos ( 



'L* + Q = fi. 



4 The conversion from black hole spin to radiative efficiency is still debated, and depends on whether or not magnetic stresses exert a 
significant torque at the last stable orbit (e.g. Krolik, Hawley, & Hirose 2004); here we adopt the standard definition for circular equa- 
torial orbits around a Kerr hole. For progradc orbits, the mass-to-encrgy conversion efficiency e = 1 — E can then be expressed as 

a = 2/(3a/3) [2- v /2/(2e - e 2 ) - (1 - e)/(2e - e 2 ) . 
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The smaller hole plunges into the large one carrying 
along its orbital constants at the LSO 5 ; given the inclina- 
tion angle, one can determine the radius of the LSO, L z , 
and E by numerically solving the system R = dR/dr = 
d 2 R/dr 2 = 0, where R is the "effective potential" gov- 
erning radial motion. The constants are bounded by their 
values for prograde (fi = 1) and retrograde (/j, = — 1) equa- 
torial orbits, 



r = r LSO /mi 



E\ 



3 + Z 2 t V(3 
E-Lso/m 2 = 



Z 1 )(3 + Z 1 + 2Z 2 ), 



1 



1/2 



3f J 



L = L L so/(m 2 mi) = ±— -=[1 + 2(3f - 2) 1/2 ] 



3V3 



= 0, 



(2) 
(3) 

(4) 
(5) 



where Z\ and Z 2 are functions of a\ only (Bardeen, Press, 
& Teukolsky 1972), and the upper (lower) sign refers to 
prograde (retrograde) orbits. According to HB, the nu- 
merical results are remarkably well fitted by 



m - 16, 



+ 2^+ 1 )(£p™ 



l&etl), 



(6) 



where £ stands for f, E, or L, and the "ret" and "pro" 
subscripts correspond to retrograde and prograde orbits, 
respectively. After capturing the smaller hole, the remnant 
MBH has mass and spin given by 



m' = toi[1 + qE(ai,fj,)] , 
S' z = m\\a x +qL z (a 1 ,^)} , 
S' ± = qmlL z (a 1 ,iJ,)y / fi- 2 - 1 . 



(7) 
(8) 
(9) 



The remnant hole is inclined at an angle A9 relative to 
the original hole, and has final spin parameter a': 



A(9 = arccos S' Z /JS' 2 



S' 2 



= JS' 2 + S'l/m' 2 



(10) 
(11) 



hole, and gas is dumped directly down the hole from the 
LSO, then the accretion of a rest mass dmo will lead to a 
change dm = E dm in the total gravitational mass and 
dS z = L z mdmo in the total angular momentum of the 
hole. The evolution of the hole spin is governed by the 
differential equation 



da 
din to 
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-4- - 2d. 



(12) 



This was integrated by Bardeen (1970) to obtain the evo- 
lution law 



TOl 



/3to 2 
\ m' 2 



1/2' 



for !!L < f W 

TOl 



for HL > r^. 

TOl 



(13) 



A hole that is initially non-rotating (f = 6) gets spun up to 
a maximally-rotating state (a' = 1) after a modest amount 
of accretion, to' = ra\\f%. A maximally-rotating hole 
(a = 1) gets spun down by retrograde accretion (f = 9) to 
a' = after to' = miy / J/2. A 180° flip of the spin of an 
extreme-Kerr hole (a = a' = 1) will occur after to' = 3toi. 

What is the spin change induced by an accretion disk 
that is inclined relative to the equatorial plane of the hole? 
Irrespective of the infalling material's original angular mo- 
mentum vector, Lense-Thirring precession will impose ax- 
isymmetry close in, with the gas accreting on either pro- 
grade or retrograde equatorial orbits. Further out, the 
disk will be warped where the transition from aligned to 
misaligned flow occurs (Bardeen & Petterson 1975). This 
suggests that accretion of material with randomly oriented 
angular momentum vectors would lead to spin-down of a 
MBH, given the larger LSO of retrograde orbits (Moderski 
& Sikora 1996). 6 However, while the material falling into 
the hole is aligned (or anti- aligned) with it, the torque 
that aligns the inner disk with the hole must ultimately 
realign the hole with the outer accretion disk (Rees 1978; 
Scheuer & Feiler 1996), thus leading to accretion via pro- 
grade equatorial orbits. The timescale for this alignment 
depends on how a warped accretion disk communicates the 
warp and evolves in time, and is given roughly by (Scheuer 
& Feiler 1996) 
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(14) 



Figures 2 and 3 show how, following a single coalescence, 
the spin of the remnant hole a! depends on the mass ratio 
and inclination cosine for two representative values of the 
initial spin oi. There is little change in spin for q < 0.05, 
while for larger mass ratios the angular momentum L at 
the LSO overwhelms the initial spin S of the hole. For 
randomly-distributed inclination angles, MBHs that are 
initially rotating with d\ = 0.6 are spun up even further 
for q > 0.125, and are spun down otherwise. Rapidly ro- 
tating holes, with cii = 0.9, are typically spun down except 
for large mass ratios and close-to-prograde orbits. 

4. SPIN CHANGES BY ACCRETION 

Accretion from a disk of gas orbiting the hole will 
cause the spin to evolve in both magnitude and orien- 
tation. If the disk is in the equatorial plane of a Kerr 

5 We have assumed here that the emission of gravitational waves during the plunge phase does not affect the binary energy and angular 
momentum. This may not be true for nearly equal mass mergers, for which there might not be a well-defined plunge phase (e.g. Pfciffcr, 
Teukolsky, & Cook 2000). 

6 This argument is similar to that given by HB for why an ensemble of binary coalescences leads to spin-down. 



where M is the accretion rate onto the hole, R$ is the 
Schwarzschild radius, R w is the location of the warp, and 
v\ and v 2 are the viscosities associated with dissipating 
motion in and out of the plane of the disk, respectively. 
For a thin disk of vertical height H, sound speed c s , and 
angular velocity O, z^/^i = (2a 2 ) -1 for 1 ^> a > H/R 7 
where v\ is given by the usual Shakura-Sunyaev prescrip- 
tion v\ = ac s H = aH 2 Q (Papaloizou & Pringle 1983). 
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MBH spin distribution 



The location of the warp in the disk (R w ) can be esti- 
mated by equating the diffusion time of the disk's warp to 
the timescale for Lense-Thirring precession, which implies 
(Natarajan & Pringle 1998) 



R w 



V2c 
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2/3 



Substituting equation (15) into equation (14) yields 



iaiign ~ 3 (~j^J (V2a 2 a 
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2/3 



2/3 



(15) 



(16) 



the equilibrium spin by transporting angular momentum 
outward. Fully relativistic magnetohydrodynamic simula- 
tions for a series of thick accretion disk models show that 
spin equilibrium is reached at a' w 0.93 (Gammic, Shapiro, 
& McKinney 2004). For simplicity, we will assume in the 
following that a MBH may be spun up to a maximum 
equilibrium value of a' = 0.998. 

5. MBH SPIN DISTRIBUTION AND EVOLUTION 

The combination of a halo merger tree and our semi- 
analytical scheme to treat the growth of MBHs and their 
dynamics is a powerful tool for tracking the evolution of 
MBH spins with cosmic time. In our merger tree, MBHs 
that undergo an accretion episode typically increase their 
mass by about one e-folding. This is required in order to 
account for the local mass density of MBHs with growth 
solely during major mergers. The distribution of frac- 
tional changes in hole mass from gas accretion is shown 
in Figure 4 for different redshift intervals. Note how, at 
all epochs, a significant fraction of accretion events leads 
to Am /mi > 2 — 3: these individual episodes will produce 
rapidly-rotating holes independent of the initial spin. 

To bracket the uncertainties and explore various sce- 
narios we have run different sets of Monte Carlo realiza- 
tions. Our "fiducial" model assumes seed holes are born 
with an initial spin of a = 0.6 (Fryer, Wooslcy, & Heger 
2002). The evolution (both in magnitude and orientation) 
of MBH spins is driven by gas accretion and black hole 
binary coalescences. Whenever a binary forms, the angle 
of inclination of the smaller hole relative to the equatorial 
plane is chosen randomly from an isotropic distribution. 
We assume that during a major merger the gas accretes 
via a thin disk; efficient alignment between the hole and 
the outer disk implies that most of the mass is accreted in 
prograde equatorial orbits (§4). 

Figure 5 shows the ensuing spin distribution of MBHs 
at different epochs. The left panel depicts the effect of 
binary coalescences alone on the spin distribution. At 
late epochs binaries with small mass ratios are common 
(Fig. 1), as some less massive holes that were wandering 
in galaxy halos at high redshift - with an orbital decay 
timescale comparable to the Hubble time - finally find 
their way to galaxy centers. In addition, our assumption 
that, following a major merger between halos, a new ac- 
cretion episode is triggered only on the more massive hole, 
probably somewhat underpredicts the mass ratio of coa- 
lescing holes. A more realistic model in which both holes 
accrete would lead to somewhat larger mass ratios, though 
this effect is likely to be mild given that the mass accreted 
per merger is comparable to the mass of the pre-existing 
hole (Fig. 4). 

Figure 5 shows that repeated captures of smaller holes 
cause a spread both towards a > 0.6 (spin up) and a < 0.6 
(spin down) in N(a). The relatively flat distribution of 
N(q) for captures in hierarchical models (Fig. 1) im- 
plies that neither spin-up nor spin-down is particularly 
favored. The holes random walk around the initial "seed" 
value, and the spin distribution retains memory of the ini- 
tial rotation black holes are endowed with at birth. For 
comparison, we also show a case in which seed holes are 

7 If the mass supply rate is supcr-Eddington, most of the inflowing gas is likely driven away in a radiation-pressure driven outflow, limiting 
the accretion rate onto the hole to ~ the Eddington rate (e.g., Shakura & Sunyaev 1973; Blandford & Begelman 2004). Thus accretion via a 
geometrically thick disk may still be compatible with the mean accretion efficiency of ~ 10% inferred by, e.g., Yu & Tremainc (2002). 



Equation (16) shows that, for thin accretion disks with 
H <C R, the alignment time is much less than the timescale 
for the mass of the hole to increase via accretion (using 
a-disk models, Natarajan & Pringle 1998 and Natarajan 
& Armitage 1999 estimate f a ii g n ~ 10 5 — 10 6 yrs for ac- 
cretion near the Eddington limit). As a result, the hole 
will align itself with the outer disk before it accretes much 
mass. This implies that most of the mass accreted by the 
hole will act to spin it up (i.e., the magnitude of the spin 
increases in time), even if the direction of the spin axis 
changes in time. This is true unless the angular momen- 
tum of the inflowing material changes significantly on a 
timescale <C t a ii g n- 

If the mass supply rate to a MBH is super-Eddington, 
accretion proceeds via a radiation pressure supported thick 
disk with H <~ R (e.g., Begelman & Meier 1982). Equa- 
tions (15) and (16) show that in this case R w ~ Rs and 
alignment is relatively inefficient, occurring on a timescale 
comparable to the timescale for changes in the hole's mass 
(the Salpeter time). It is possible that much of the growth 
of MBHs at high redshift occurs via a radiation pressure 
supported thick disk, rather than a thin accretion disk. 
In fact, it would be surprising if the mass supply rate were 
precisely ~ Eddington (required for a thin disk) during the 
entire growth of MBHs. It is perhaps more likely that the 
mass supply rate is initially larger in the dense gas-rich 
environments of high redshift galaxies and then decreases 
with time as the galaxy is assembled. 

In the extreme case in which alignment is always rela- 
tively inefficient, as would be expected for accretion via a 
thick disk, we can generalize equations (7), (8), and (9) 
to determine the angular momentum evolution under gas 
accretion as follows: 



dm = dm E(ai, fi) , 

dS z =midm L z (ai,jj l ) , 

dS± = midmoL z (ai, fj,) \J /z -2 — 1 



(17) 
(18) 
(19) 



This equation can be integrated numerically to yield the 
final orientation and magnitude of the hole spin after each 
accretion episode, for a given inclination. 

Finally, we note that Thorne (1974) showed that the 
radiation emitted by the disk and swallowed by the hole 
produces a counteracting torque, which prevents spin up 
beyond a! — 0.998. Magnetic fields connecting material 
in the disk and the plunging region may further reduce 
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born non-spinning instead (dashed histogram in the left 
panel). Again, at all redshifts, the spin distribution re- 
mains peaked around a = 0. The right panel in Figure 
5 includes spin changes both by captures and by gas ac- 
cretion. Gas accretion dominates the spin evolution over 
coalescences (the initial mass in seed holes is a tiny frac- 
tion of the mass density observed today). At all epochs, 
nearly all MBHs are spun-up by accretion to a w 1. 

Figure 6 shows the spin evolution history of two MBHs 
that end up in massive halos with M/j = 10 12 M Q at z = 0. 
Note how gas accretion (lower panel) efficiently spins the 
holes up, while binary coalescences lead to both spin-up 
and spin-down. Figures 7 and 8 explore the uncertainties 
introduced by some of our assumptions. Figure 7 shows 
the evolution of the spin under binary coalescences alone, 
comparing two examples where we placed seed holes in 
3<T and 3.5c peaks at high rcdshift. The former exam- 
ple has a total of <~ 50 times more seed holes, leading to 
<~ 10 times more binary coalescences; 8 the total mass in 
seed holes is now a larger fraction of the mass density in 
MBHs observed today. Despite the significant increase in 
the number of binary coalescences, the spin evolution is 
reasonably similar in the two examples. The most notable 
difference is that there are fewer rapidly rotating holes at 
low redshift in the case with more seed holes, because of 
spin-down by capture of smaller companions at late times 
(see Fig. 1). The overall effect is, however, rather mild, 
and the spin distribution still retains significant memory 
of the initial conditions; binary coalescences alone thus do 
not lead to a significant systematic spin-up or spin-down 
of MBHs. 

Figure 8 shows an example where alignment of the hole 
and the outer disk is inefficient, as would be expected if 
accretion is via a geometrically thick disk (§4). We as- 
sume that the initial orientation between the black hole's 
spin and the accretion disk rotation axis is random, and 
integrate equations (17)-(19) for every accretion episode, 
i.e. we treat binary coalescences and gas accretion in a 
similar way. This would seem to favor spin-down of the 
MBHs, and yet even in this case most of the holes are 
rotating quite rapidly at all epochs. The reason is that 
the accreted mass is typically larger than the hole's mass 
(Fig. 4). Most individual accretion episodes thus produce 
rapidly-rotating holes independent of the initial spin. 

6. SUMMARY AND DISCUSSION 

We have computed the expected distribution of MBH 
spins and its evolution with cosmic time in the context of 
hierarchical structure formation theories. A subset of the 
current authors have previously described models for the 
birth, growth, and dynamics of MBHs that traces their 
build-up far up the dark halo "merger tree" in a ACDM 
cosmology (see, e.g., Volonteri et al. 2003). Here we have 
extended this work to follow the combined effects of black 
hole-black hole coalescences and accretion from a gaseous 
disk on the magnitude and orientation of MBH spins. 

We find that binary coalescences cause no significant 
systematic spin-up or spin-down of MBHs: because of the 
relatively flat distribution of MBH binary mass ratios in hi- 
erarchical models (Fig. 1), the holes random- walk around 

8 The nonlinear increase in the number of binary coalescences in going from 3.5 to 3 a peaks is due to several factors: (1) there are a larger 
number of halo minor mergers that never form MBH binaries because of the long dynamical friction timescales, and (2) a larger number of 
holes arc ejected from their host galaxies by triple interactions and gravitational recoil. 



the spin parameter they are endowed with at birth, and 
the N(a) distribution retains significant memory of the 
initial rotation of "seed" holes (Figs. 5 and 7). 

In our models, accretion, not binary coalescences, dom- 
inates the spin evolution of MBHs. Accretion can lead to 
efficient spin-up of MBHs even if the angular momentum 
of the inflowing material varies in time. This is because, 
for a thin accretion disk, the hole is aligned with the outer 
disk on a timescale that is much shorter than the Salpeter 
time (eq. 16; Natarajan & Pringle 1998), leading to ac- 
cretion via prograde equatorial orbits. As a result, most 
of the mass accreted by the hole acts to spin it up, even if 
the orientation of the spin axis changes in time. For a geo- 
metrically thick disk, alignment of the hole with the outer 
disk is much less efficient, occurring on a timescale com- 
parable to the Salpeter time. In this case we still find that 
most holes are rotating rapidly (Fig. 8). This is because, 
in any model in which MBH growth is triggered by major 
mergers, every accretion episode must typically increase a 
hole's mass by about one e-folding to account for the local 
MBH mass density and the tobh — relation. Most in- 
dividual accretion episodes thus produce rapidly-rotating 
holes independent of the initial spin. 

Under the combined effects of accretion and binary co- 
alescences, we find that the spin distribution is heavily 
skewed towards fast-rotating Kerr holes, is already in place 
at early epochs, and does not change significantly below 
rcdshift 5. As shown in Figure 9, about 70% of all MBHs 
are maximally rotating and have mass-to-energy conver- 
sion efficiencies approaching 30%. Note that if the equi- 
librium spin attained by accreting MBHs is lower than the 
value of a — 0.998 used here, as in the thick disk MHD 
simulations of Gammic ct al. (2004) where a w 0.93, then 
the accretion efficiency will be lower as well, « 17 % for 
a w 0.93. Even in the conservative case where accretion 
is via a geometrically thick disk (and hence the spin/disk 
alignment is relatively inefficient) and the initial orienta- 
tion between the hole's spin and the disk rotation axis is 
assumed to be random, we find that most MBHs rotate 
rapidly with spin parameters a > 0.8 and accretion effi- 
ciencies e >12%. As recently shown by Yu & Tremaine 
(2002), Elvis, Risaliti, & Zamorani (2002), and Marconi 
et al. (2004), a direct comparison between the local MBH 
mass density and the mass density accreted by luminous 
quasars shows that quasars have a mass-to-energy conver- 
sion efficiency e > 0.1 (a simple and elegant argument 
originally provided by Soltan 1982). This high average ac- 
cretion efficiency may suggest rapidly rotating Kerr holes, 
in agreement with our findings. 

In our models, there is a weak trend for the most mas- 
sive and the least massive holes to have slightly lower spin 
parameters (efficiencies), the former because they experi- 
ence more binary coalescences during their lifetime, and 
the latter because they experience a smaller number of 
accretion episodes. 

One way to avoid rapid rotation and produce slowly ro- 
tating holes is to assume "chaotic feeding" in which small 
amounts of material, with Am <C tobh, arc swallowed 
by the hole in successive accretion episodes with random 
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orientations (e.g., Moderski and Sikora 1996). The con- 
straints on such a model appear, however, to be strin- 
gent. Specifically, for a thin disk the angular momen- 
tum of the inflowing material must vary on a timescale 
fti iaiign <C ^Saipetcr- Otherwise alignment of the hole with 
the outer disk ensures that most of the mass accreted by 
the hole spins it up. The required timescale for angular 
momentum reversals is much shorter than that expected 
in merger-initiated nuclear activity, but could in princi- 
ple be provided by (say) accretion of molecular clouds on 
random orbits. 

In this context, an issue which must be investigated fur- 
ther is the alignment of the BH spin axis and the angular 
momentum axis of the outer disk when the initial angle 
is more than 90 degrees. In this case the inner disk ends 
up "anti-aligned" instead of "aligned"; to fully align the 
hole with the outer disk requires the inner disk to "flip" 
by 180 degrees. If the timescale for alignment is as short 
as we have assumed for the thin disk geometry, then the 
initial period of "anti-alignment" is probably not impor- 
tant, but if the timescale is of order the Salpeter timescale 



successive accretion events could cancel out. 

To explain the dichotomy between radio-quiet and 
radio-loud quasars, Wilson & Colbert (1995) suggested 
that in all radio-quiet quasars the MBH is slowly rotating. 
Since most holes rotate rapidly in our models, our results 
suggest instead that the spin of a MBH is not a necessary 
and sufficient condition for producing a radio-loud quasar. 
Independent evidence for this conclusion comes from ob- 
servations of X-ray binaries, where jets are observed to be 
present only in the hard/low X-ray spectral state, but not 
in the soft/high state (Fender 2001). This favors a model 
in which the mode of accretion, perhaps together with the 
spin, determines the radio loudness of an accreting black 
hole (e.g., Rees et al. 1982; Meier 2001). 
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Fig. 1. — Normalized distribution of mass ratios of coalescing MBH binaries at three different epochs. Note that at low redshift MBHs 
typically capture much smaller companions. 




Fig. 2. — Final spin a! of the remnant MBH as a function of inclination angle fi = cos$, for different binary mass ratios q = m2/m\. Before 
coalescence, the larger hole has spin &i = 0.6. The remnant is typically spun down for q < 0.125, and is spun up otherwise. 
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Fig. 3. — Same as Fig. 2, but with ai = 0.9. For most of the parameter space, the remnant MBH is spun down. 
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Fig. 4. — Distribution of fractional changes in the mass of a MBH as a result of the accretion of material after a major halo merger. 
Different redshift intervals arc shown. 
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Fig. 5. — Distribution of MBH spins in different redshift intervals. Left panel: effect of black hole binary coalescences only. Solid histogram: 
"seed" holes are born with a = 0.6. Dashed histogram: "seed" holes are born non-spinning. Right panel: spin distribution from binary 
coalescences and gas accretion. "Seed" holes are born with a = 0.6, and are efficiently spun up by accretion via a thin disk. 
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Fig. 6. — The spin evolution of two MBHs, "A" (left panels) and "B" (right panels), that end up in massive halos with = 10 12 Mq at 
z = 0. The initial spin of "seed" holes is a = 0.6. Upper panel: the spin of black holes is modified by binary coalescences only. Lower panel: 
the spin evolution is driven by black hole binary coalescences and gas accretion. The dots mark the epoch when two MBHs coalesce. Note 
the intermittent times when halos are devoid of a central MBH, after an ejection due to radiation recoil. The halo (re)gains a central MBH 
after the original one falls back in or the halo experiences a merger with another halo, whose MBH becomes the central hole of the newly 
formed galactic system. 
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Fig. 7. — The spin distribution for MBHs that end up in massive halos with = 10 12 Mq at z = 0. Spins evolve under by binary 
coalescences alone. Left panel: "seed" holes are placed in 3cr peaks at z = 20. Right panel: "seed" holes are placed in 3.5<r peaks at z = 20. 
In both cases black holes are born with a = 0.6. The model in the left panel has ~ 10 times more coalescences than that in the right panel, 
but the resulting spin distributions are quite similar. 
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Fig. 8. — Distribution of MBH spins in different rcdshift intervals. The spin is modified by binary coalescences and by accretion; alignment 
of the hole and the outer disk is assumed to be inefficient, as would be expected for accretion via a geometrically thick disk. The initial 
orientation between hole's spin and disk angular momentum is assumed to be random. 
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Fig. 9. — Distribution of accretion efficiencies, e = 1 — E, in different redshift intervals (assuming that the energy radiated is the binding 
energy at the LSO). The spin distribution from binary coalescences and gas accretion has been calculated assuming the holes accrete via a 
thin disk on progradc equatorial orbits. 



